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Abstract.  Recommendation  VII  of  Resolution  A4  of  the  XXIst  General  Assembly  of  the  International  Astronomical  Union 
(IAU  1992)  states,  in  part,  “...  that  the  principal  plane  of  the  new  conventional  celestial  reference  system  be  as  near  as  possible 
to  the  mean  equator  at  J2000.0  and  that  the  origin  in  this  principal  plane  be  as  near  as  possible  to  the  dynamical  equinox 
of  J2000.0,  ...”  The  resulting  International  Celestial  Reference  System  (ICRS),  however,  has  a  small,  but  significant,  offset 
requiring  a  rotation  matrix.  The  solutions  for  the  offset  between  the  mean  dynamical  pole  of  the  Earth  at  J2000.0  and  the  pole 
of  the  ICRS  determined  by  Lunar  Laser  Ranging  (LLR)  and  Very  Long  Baseline  Interferometry  (VLBI)  differ  by  several  <x. 
Similarly,  two  different  definitions  have  traditionally  been  used  for  the  position  of  the  mean  equinox.  Which  of  these  poles  and 
equinoxes  should  be  used  is  application  dependent.  We  have  shown  how  the  rotation  matrix  for  the  rotation  from  the  mean 
dynamical  equator  and  equinox  at  J2000.0  to  the  ICRS  changes  depending  on  the  various  assumptions  made  in  constructing  it. 
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1.  Introduction 

Resolution  A4  VII  of  the  XXIst  General  Assembly  of  the  International  Astronomical  Union  (IAU  1992)  recommended  the  adop¬ 
tion  of  a  quasi-inertial  reference  system  based  on  the  Very  Long  Baseline  Interferometry  (VLBI)  observations  of  extragalactic 
radio  sources.  This  reference  system  was  to  be,  “...  computed  initially  for  the  equator  and  equinox  J2000.0  using  the  best  avail¬ 
able  values  of  the  celestial  pole  offset  with  respect  to  the  IAU  expressions  for  precession  and  nutation  ...”  The  final  reference 
system,  the  ICRS,  was  adopted  in  1997  in  resolution  B2  of  the  XXIIIrd  assembly  of  the  IAU  (IAU  1997).  Since  that  time,  im¬ 
proved  measurements  have  shown  that  there  is  a  small,  but  significant,  difference  between  the  orientation  of  the  ICRS  and  that 
of  the  best  estimate  for  the  mean  dynamical  equator  and  equinox  of  J2000.0  (henceforth  J2000.0).  The  difference  between  the 
orientations  of  the  ICRS  and  J2000.0  amounts  to  a  few  tens  of  milliarc seconds  (mas).  This  difference  can  be  significant  for  those 
working  at  the  highest  levels  of  precision.  For  example,  application  of  precession  and  nutation  to  the  ICRS  rather  than  J2000.0 
leads  to  additional  errors  that  grow  to  the  same  order  of  magnitude  as  the  difference  in  orientation  between  the  ICRS  and  J2000.0 
(Wallace,  priv.  comm.  2002).  We  have  determined  the  rotation  matrix  from  J2000.0  to  the  ICRS  using  the  current  best  estimates 
for  these  two  reference  systems. 

2.  Derivation  of  the  rotation  matrix 

The  coordinate  system  of  J2000.0  can  be  described  using  a  set  of  Cartesian  coordinate  axes,  Xj2ooo,  Yj2ooo-  Znooo,  where  the 
positive  Xjoooo  axis  points  towards  the  equinox,  the  positive  Yqooo  axis  points  towards  the  direction  in  the  equator  of  the  system 
that  is  90°  to  the  east  of  the  positive  Xj2ooo  axis,  and  the  positive  Zj2ooo  axis  points  in  the  direction  Xj2ooo  ><  Yj2ooo-  Similarly, 
the  ICRS  coordinate  system  can  be  described  using  a  set  of  Cartesian  coordinate  axes,  Xicrs,  Yicrs,  Zicrs,  where  the  positive 
Xicrs  axis  points  towards  the  origin  of  the  ICRS  coordinate  system,  the  positive  Yicrs  axis  points  towards  the  direction  in 
the  equator  of  the  system  that  is  90°  to  the  east  of  the  positive  Xicrs  axis,  and  the  positive  Zicrs  axis  points  in  the  direction 
Xicrs  x  Yicrs- 

The  transformation  of  a  vector  x  from  the  J2000.0  coordinate  system,  xj2ooo>  to  the  ICRS,  Xicrs,  by  a  rotation  matrix  R  is: 
xicrs  =  R  xnooo  =  R>p  Re  R<p  XJ2000  (1) 
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J.  L.  Hilton  and  C.  Y.  Hohenkerk:  Rotation  matrix  from  J2000.0  to  the  ICRS 


Fig.  1.  The  position  of  the  pole  of  J2000.0,  Pkooo.o.  with  respect  to  the  pole  of 
the  ICRS,  Picrs-  The  relative  position  is  given  by  the  position  angle,  PA  and  the  distance  d, 
or,  alternatively  by  the  Cartesian  coordinates  ex  and  ey  where  Xp  and  Yp  are  the  coordinate 
axes  described  by  Feissel  &  Mignard  (1998). 


where  R^,  Rg ,  and  R,/,  are  the  rotation  matrices  for  the  three  Eulerian  rotations.  The  Euler  angles  are:  <//,  the  rotation  from  the 
node  of  the  two  systems  to  the  origin  of  the  second  system  around  the  polar  axis  of  the  second  system;  0,  the  inclination  of  the 
second  system  with  respect  to  the  first  system;  and  <p ,  the  initial  rotation  around  the  polar  axis  of  the  first  system  from  its  origin 
to  the  node  of  the  two  systems.  The  individual  Euler  rotation  matrices  are: 


Rip 

Re 

R<p 


v 

/ 


cos  ip 
-  sin  0- 


0 


sin  0- 

COSl/r 

0 


O' 

0 


10  O' 
0  cos  6  sin  6 
0  -  sin  8  cos  6  / 

cos  (f>  sin  (j)  0  ' 
-  sin  (A  cos  (j)  0 


and 


0  0  1  J 


(2) 


The  difference  in  position  of  the  J2000.0  pole  with  respect  to  the  pole  of  the  ICRS  is  given  as  an  offset,  as  published  by  Feissel 
&  Mignard  (1998)  and  the  International  Earth  Rotation  Service  (IERS)  (IERS  Annual  Reports).  This  offset  uses  an  Xp  -  Yp 
coordinate  system,  where  the  positive  Xp  direction  is  the  6  hr  meridian  of  the  ICRS  and  the  positive  Yp  direction  is  the  12  hr 
meridian,  as  shown  in  Fig.  1.  Note  that  this  coordinate  system  is  rotated  by  90°  from  the  projection  of  the  ICRS  Xicrs  -  Yicrs 
plane  onto  the  Xp  -  Yp  plane.  PA  is  the  position  angle,  measured  from  the  positive  Xp  axis,  of  J2000.0  with  respect  to  the  ICRS 
and  d  is  the  angular  distance  between  the  pole  of  the  ICRS  and  the  pole  of  J2000.0.  Since  the  Xp  -  Yp  plane  is  rotated  with  respect 
to  the  projection  of  the  Xicrs  -  Yicrs  plane,  aicRS  -  PA  =  90°  where  q-icrs  is  the  ICRS  meridian  coordinate1.  Alternatively,  the 
position  of  the  pole  of  J2000.0  with  respect  to  the  ICRS  can  be  given  in  the  Xp  -  Yp  plane  by  the  Cartesian  coordinates  ex  and 
ey,  The  relations  between  the  Cartesian  and  polar  coordinates  are: 


ex  =  d  cos  PA  and  ey  =  d  sin  PA.  (3) 

The  relative  positions  of  the  origin  of  J2000.0  on  its  equator,  and  the  origin  of  the  ICRS  on  its  equator,  are  shown  in  Fig.  2.  The 
origins  of  both  the  ICRS  and  J2000.0  are  supposed  to  coincide  with  their  respective  equinoxes.  Chapront  et  al.  (2002),  however, 
determined  that  the  actual  origin  of  the  ICRS  is  offset  from  its  equinox  by  A o.  The  angle  eo  is  the  obliquity  of  the  ecliptic  on  the 
ICRS,  that  is  the  angle  between  the  ICRS  equator  and  the  mean  ecliptic  as  determined  by  Chapront  et  al.  from  LLR  observations 
(23° 26' 2 1741 100  +  07 00005).  The  angle  ytJ  is  the  separation  between  the  projection  of  the  node  of  the  mean  ecliptic  on  the  ICRS 
onto  the  equator  of  J2000.0  and  the  equinox  of  J2000.0  shown  in  Fig.  2.  This  angle  is  given  by: 


tan  yy  =  tan  At [t  cos  eo  (4) 

1  We  use  the  term  meridian  coordinate  instead  of  right  ascension  for  the  longitudinal  coordinate  of  the  ICRS  because  the  ICRS  is  insensitive 
to  the  equinox,  the  traditional  origin  of  right  ascension,  and  its  origin  is  not  aligned  with  the  node  of  its  equator  on  the  mean  ecliptic. 
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_ h. _ 

J2000.0  equator 

1  ^^yI2„0„(J2OOO.O) 

CRS  equator 

Y200o(ICRS) 

o(ICRS) 

A  0 

Fig.  2.  The  relative  positions  of  the  origins  of  the  ICRS  and  J2000.0  at  J2000.0.  The  origin  of  J2000.0  is  the  same  as  the  position  of  the  inertial 
equinox,  T-iooo^OOO.O).  The  origin  of  the  ICRS,  however,  is  offset  from  its  inertial  equinox,  y^/ICRS),  by  A o.  The  projection  of  the  offset 
between  y(000(ICRS)  and  y^ooo^OOO.O)  on  the  ICRS  equator  is  yy.  And  the  separation  between  the  two  equators  at  y2000(ICRS)  is  the  same  as 
the  y  coordinate  of  the  separation  of  the  two  poles,  ey.  The  angle  is  positive  when  it  goes  from  left  to  right  or  bottom  to  top  in  the  figure. 

Table  1.  Recent  determinations  of  the  angles  between  the  ICRS  and  J2000.0. 


Authors 

Method  of 

A«o 

Ai// 

A  oa 

observation 

(mas) 

(mas) 

(mas) 

Chapront  et  al.  (2002) 

LLR 

-5.36  ±0.10 

-44.5  ±  0.3 

55.42  ±0.11 

Mathews  et  al.  (2002) 

VLB  I 

-6.8192 

-41.775 

- 

“  The  angle  A o  is  measured  from  the  node  of  the  ICRS  with  the  mean  ecliptic  of  J2000.0  to  the  origin  of  the  ICRS.  Thus,  it  has  the  opposite 
sign  from  that  published  in  Chapront  et  al.  (2002),  which  is  measured  in  the  opposite  direction  (see  text  for  explanation). 


Table  2.  Rotation  parameters  between  the  ICRS  and  J2000.0. 


Authors 

Method  of 

e« 

1 

< 

observation 

(mas) 

(mas) 

(mas) 

(mas) 

Chapront  et  al.  (2002) 

LLR 

-5.36  ±0.54 

17.70  ±0.20 

40.83  ±  0.46 

14.6  ±  0.5 

Mathews  et  al.  (2002) 

VLB  I 

-6.819  ±0.015 

16.6171  ±0.0080 

38.328  ±0.018 

- 

where  Ai jj  is  the  angle  between  the  pole  of  the  ICRS  and  the  pole  of  J2000.0  as  seen  from  the  ecliptic  pole.  Both  yy  and  At//  are 
small  angles,  so  they  can  be  expanded  to  first  order.  The  error  in  yy  introduced  by  using  the  expansions  is  ~  1 0  1 5  radians,  much 
smaller  than  the  uncertainty  in  At/t  of  10  y-10  n  radians. 

Note  that  yy  is  defined  such  that  the  initial  rotation  about  the  polar  axis  is  from  the  node  of  the  mean  ecliptic  on  the  ICRS 
to  the  node  of  the  mean  equator  of  J2000.0,  not  from  the  origin  of  the  ICRS  to  the  node  of  the  mean  equator  of  J2000.0.  The 
result  is  that  the  values  of  A o  in  Table  1  and  ytJ  in  Table  2  have  the  opposite  signs  than  they  would  if  they  were  measured  from 
the  origin  of  the  ICRS  and  mean  equinox  of  J2000.0,  respectively. 

Using  these  angles,  the  three  Eulerian  rotation  angles  are  then: 


if>  =  PA  +  A o,  8  —  d,  and  t f/  =  -PA  -  yy. 


(5) 


Simplification  of  the  elements  of  the  rotation  matrix  by  approximating  the  sine  and  cosine  functions  is  useful  in  understanding 
the  rotation  between  the  two  coordinate  systems. 

Expanding  sin  x  and  cos  x  to  first  order  implies  that,  for  elements  in  the  rotation  matrix  that  contain  sine  of  a  small  angle, 
the  uncertainty  in  that  element  is  <x~  ~  SI  •  cr\  while  those  elements  containing  only  the  cosine  function  have  an  uncertainty  of 
cr“  ~  'Ll/1  ■  (Xx.  The  small  angles  in  the  rotation  matrix  are  ~  1 0  mas  (5  x  1 0  8  rad)  with  uncertainties  ~0. 1  mas  (5  x  10~10  rad). 
Thus,  the  uncertainty  of  those  elements  that  contain  a  sine  function  is  ~5  x  10  10 ,  and  only  need  to  be  expanded  to  first  order.  On 
the  other  hand,  those  elements  that  contain  only  cosine  functions  have  uncertainties  ~  1 0  18  and  need  to  be  expanded  to  second 
order. 

When  the  values  for  the  angles  are  substituted  into  the  rotation  matrix  and  the  simplifications  are  carried  out,  R  becomes: 


,  (A  o-yy)2 
1  2  2 

-(A o  -  yy)  1 
<% 


A O  -  Jy  ~€y 

(Ao-rj,)2  el 

2  2 


(6) 


The  appearance  of  ey  in  the  R\\  and  R)i  elements,  ex  in  the  R23  and  R22  elements  and  the  change  in  sign  between  ex  and  ey  in  the 
third  column  and  third  row  appear  to  be  counterintuitive.  The  reason  that  these  elements  appear  in  this  order  with  these  signs  is 
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Fig.  3.  The  position  of  the  poles  of  the  ICRS,  Picrs,  and  J2000.0,  P.;20oo.o 
with  respect  to  the  pole  of  the  ecliptic,  C.  The  angle  between  the  two  poles 
is  Ai/r  and  the  difference  in  their  obliquities  is  Ae0  =  6ICrs  -  £j20oo.o- 


due  to  the  rotation  of  the  Xp  -  Yp  plane  with  respect  to  the  projection  of  the  Xicrs  -  Yicrs  plane.  A  position  along  the  positive  Xp 
is  the  equivalent  of  a  position  along  the  positive  Yicrs  axis,  while  a  position  along  the  positive  Yp  axis  is  the  equivalent  of  a 
position  along  the  negative  Xicrs  axis. 

Note  also,  although  for  parameters  ex,  ey,  yy,  and  Ao  appear  in  (6)  only  three  of  them  are  independent.  The  parameter  yy  is 
related  to  Ai jj  by  (4)  and  ex  and  ey  are  related  to  At//  and  Aeo,  the  difference  in  the  obliquity  of  the  ICRS  pole  and  the  J2000.0  pole, 
via  (3)  and  (7)  below. 


3.  Parameter  values 


Mathews  et  al.  (2002)  estimate  the  offset  between  the  pole  of  J2000.0  and  the  pole  of  the  ICRS  from  VLBI  observations.  The 
ICRS  is  also  defined  by  VLBI  observations  of  extragalactic  objects  (Ma  et  al.  1998).  The  relative  positions  of  the  two  poles,  shown 
in  Fig.  3,  are  given  by  A \Jj  (4)  and  Aeo-  Mathews  et  al.  do  not  provide  errors  in  their  values;  however,  other  recent  determinations 
of  these  parameters,  such  as  Shirai  &  Fukushima  (2001),  suggest  that  the  uncertainty  is  about  0.01  mas  for  both  parameters.  This 
uncertainty  will  be  adopted  for  the  Mathews  et  al.  parameters.  From  the  spherical  triangle  in  Fig.  3,  the  relations  between  At//  and 
Aeo  and  the  angles  d  and  PA  are: 


ej2000  -  cicrs  +  Aeo 

cos  d  -  cos  eicRs  cos  enooo  +  sin  eicRs  sin  ejaooo  cos  A t/r 

.  sinAi^  . 

sin  PA  =  — — -  smeJ2ooo- 
sind 


(7) 


VLBI,  however,  is  insensitive  to  the  position  of  the  equinox  and  the  origin  of  the  ICRS.  The  best  estimate  for  these  positions  is 
determined  by  Chapront  et  al.  (2002)  from  LLR  data.  Table  1  gives  the  values  for  the  offsets  between  J2000.0  and  ICRS  found 
by  Chapront  et  al.  (LLR)  and  Mathews  et  al.  (VLBI). 

There  are  at  least  two  possible  ways  of  incorporating  the  equinox  offset  into  the  rotation  matrix.  The  method  adopted  by 
the  IAU  SOFA  software  and  the  IERS  (Wallace  2002;  Capitaine  et  al.  2003),  is  to  use  the  offset,  Ao,  of  the  origins  purely  as 
determined  by  LLR  and  the  position  of  the  pole  as  determined  by  VLBI  observations.  The  rotation  matrix  is: 


Rihrs  - 


0.99999999999999332  +  4  x  lO"17  (7.2  ±  0.2)  x  10~8  -(8.056  ±  0.008)  x  10~8 

-(7.2  ±  0.2)  x  10~8  0.99999999999999332  +  4  x  10^17  -(3.306  ±  0.004)  x  10-8 

(8.056  ±  0.008)  x  10~8  (3.306  ±  0.004)  x  IQ-8  0.999999999999996208  +  4  x  10~18  / 


(8) 


Aside  from  the  possible  systematic  error  caused  by  the  misalignment  of  the  origin,  the  accuracy  in  the  rotated  position  varies 
from  about  0.5  mas  near  the  equator  to  about  0.02  mas  near  the  poles  The  variation  in  the  accuracy  is  caused  by  the  much  lower 
accuracy  of  the  separation  between  origin  of  the  ICRS  and  the  mean  equinox  of  J2000.0  when  compared  to  the  accuracy  of  the 
separation  between  the  poles. 

An  alternative  method  is  to  use  the  value  of  A o  found  from  the  LLR  observations  with  the  VLBI  value  for  yy  to  give  Ao  -  yy. 
This  method  assumes  that  the  offset  of  the  ICRS  origin  from  the  ICRS  equinox  is  the  same  in  both  methods.  VLBI  observations, 
like  LLR  observations,  take  place  in  an  inertial  framework;  thus,  the  value  of  ytJ  determined  from  VLBI  observations  assumes 
an  inertial  mean  equinox.  The  value  of  A o  -  yy  is  then  (55.42  ±  0.11)  —  (38.33  +  0.02)  =  17.09  ±0.11  mas  in  comparison  to  the 
Chapront  et  al.  value  of  14.6  +  0.5  mas  used  in  deriving  (8).  Thus, 


Rn  =  - R2l  =  (8.29  ±  0.05)  x  10~8. 


(9) 


The  difference  between  the  Mathews  et  al.  VLBI-based  parameters  and  the  Chapront  et  al.  LLR-based  parameters  means  that 
there  is  a  possible  inconsistency  in  the  R n  and  AS|  elements  of  the  rotation  matrix  that  could  leave  a  misalignment  in  the  rotated 
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origin  of  up  to  3  mas.  Thus,  the  user  also  has  to  be  aware  that  there  may  be  a  misalignment  of  up  to  3  mas  in  the  rotated  origin, 
because  the  method  of  determining  the  position  of  the  origin  is  not  wholly  consistent  with  that  of  determining  the  position  of  the 
pole. 

On  the  other  hand,  the  node  of  the  mean  ecliptic  in  the  ICRS,  yL00(ICRS),  and  the  equinox,  Tiooo^OOO.O),  shown  in  Fig.  2 
are  specifically  for  the  inertial  definition  of  the  ecliptic  as  described  in  Standish  (1981).  The  inertial  definition  of  the  mean 
equinox  is  the  one  that  appears  in  the  IERS  Conventions  for  2000  (IERS,  2003).  However,  older  ephemerides,  such  as  DE200  are 
aligned  to  the  rotating  definition  of  the  mean  equinox  Folkner  et  al  (1994).  Standish  (1981)  found  a  difference  in  the  position  of 
the  rotating  minus  inertial  ecliptic  of  -93.66  mas  at  J2000.0.  Both  the  78+10  mas  counterclockwise  offset  of  the  rotating  equinox 
of  the  DE200  ephemerides  from  the  origin  of  the  ICRS  found  by  Folkner  et  al.,  and  the  14.6  ±  0.5  mas  clockwise  offset  of  the 
inertial  equinox  found  by  Chapront  et  al.  are  consistent  with  Standish’s  estimate.  Although  Standish  also  found  a  difference  in 
the  rotating  minus  inertial  obliquity  of  3.34  mas  at  J2000.0,  this  does  not  affect  the  relative  position  of  the  pole  of  J2000.0  with 
respect  to  the  pole  of  the  ICRS,  as  the  difference  in  obliquity  is  caused  solely  by  a  change  in  the  position  of  the  pole  of  the  mean 
ecliptic.  Thus,  the  sole  effect  of  using  the  rotating  definition  of  the  mean  equinox  rather  than  the  inertial  definition  is  to  change 
the  value  of  A o  to 


A o  =  55.42  -  93.66  =  -38.24  mas. 


(10) 


Thus,  the  R\2  and  R2\  elements  to  become: 

Rv  =  -R2l  =  -(38.3  +  0.2)  x  1(T8  (11) 

if  both  Ao  and  yy  are  determined  by  LLR.  If  A o  is  determined  by  LLR  and  yy  is  determined  from  VLBI  then 

Rv  =  -R2l  =  -(37.1  ±0.2)  x  10-8.  (12) 


The  source  of  this  difference  is  the  same  possible  3  mas  misalignment  of  the  equinox  caused  by  the  difference  in  the  VLBI  and 
LLR  solutions  for  yy  that  causes  the  difference  between  (9)  and  the  R\2  and  R2\  elements  in  (8).  Note  also  that  the  signs  of 
the  R\2  and  R  \  elements  change  indicating  that  the  rotating  equinox  is  on  the  opposite  side  of  the  origin  of  the  ICRS  from  the 
inertial  equinox. 

The  uncertainty  in  the  position  of  the  rotating  mean  equinox  with  respect  to  the  inertial  mean  equinox  is  assumed  to  be  the 
same  as  the  uncertainty  of  the  frame  tie  between  DE405  and  the  ICRS  of  1  mas  (Standish  1998).  Thus,  the  accuracy  of  the 
rotation  using  the  rotating  definition  of  the  mean  equinox  is  approximately  1.2  mas. 

Chapront  et  al.  (2002)  also  determined  values  for  At//  and  Aeo,  but  the  uncertainty  in  their  values  is  about  an  order  of  magnitude 
greater  than  those  found  from  VLBI.  Note  that  in  Table  1  the  orientation  of  the  J2000.0  pole  with  respect  to  the  ICRS  pole  found 
by  LLR  is  significantly  different  than  that  found  by  VLBI.  The  rotation  parameters  derived  from  Chapront  et  al.  and  Mathews 
et  al.  are  given  in  Table  2. 

The  rotation  matrix  from  J2000.0  to  the  ICRS  based  on  the  Chapront  et  al.  (LLR)  parameters  is: 


(  0.9999999999999938  ±  2  x  KT16 


^LLR  = 


-(7.1  ±0.2)  x  1(T8 
(8.6  ±  0.3)  x  1(L8 


(7.1  ±0.2)  x  1(T8 
0.9999999999999938  ±  2  x  1CT16 
(2.6  ±  0.1)  x  1(T8 


-(8.6  ±  0.3)  x  1(L8 
-(2.6  ±  0.1)  x  1(T8 
0.99999999999999598+  5  x  1(T17, 


(13) 


This  rotation  matrix  has  the  advantage  of  being  determined  from  a  wholly  consistent  set  of  parameters.  Its  accuracy  is  limited  by 
the  accuracy  of  the  LLR  observations.  Thus,  the  resulting  position  after  rotation  should  be  accurate  to  approximately  0.9  mas. 
Table  3  presents  the  angles  ex,  ey,  yy,  and  A o  for  use  with  (6),  for  determining  each  of  the  five  rotation  matrices  above. 


4.  Conclusions 

The  orientation  of  ICRS,  as  stipulated  in  Recommendation  VII  of  Resolution  A4  of  the  XXI  General  Assembly  of  the 
International  Astronomical  Union,  was  to  be  aligned  as  closely  as  possible  with  the  mean  celestial  ephemeris  pole  and  equinox 
of  J2000.0.  Since  its  inception,  however,  better  observations  have  shown  a  small  but  significant  offset  between  the  ICRS  and  the 
mean  equator  and  equinox  of  J2000.0.  The  parameters  that  define  this  offset  have  been  determined  by  Chapront  et  al.  (2002) 
using  LLR  observations,  and  Mathews  et  al.  (2002)  from  VLBI  data.  These  two  sets  of  parameters  are  significantly  different,  and 
the  VLBI  data  are  insensitive  to  the  position  of  the  equinox  and  the  origin  of  the  ICRS. 

There  are  also  two  different  definitions  for  the  mean  equinox  of  J2000.0  (Standish  1981),  one  defined  using  an  inertial  ecliptic 
and  the  other  based  on  the  rotating  ecliptic.  The  IAU  (1991 )  recommended  the  adoption  of  the  rotating  definition,  while  the  IERS 
(2003)  used  the  inertial  definition  in  setting  up  the  ICRS. 

We  have  thus  determined  five  scenarios  to  transform  positions  from  J2000.0  to  the  ICRS.  The  parameters  of  which  are 
summarized  below  and  listed  in  Table  3. 

The  matrix  /?iers,  (8),  uses  an  offset  of  the  origin  of  the  ICRS  from  the  mean  inertial  equinox  of  J2000.0  determined  wholly 
from  LLR  data  and  an  offset  for  the  poles  determined  wholly  from  VLBI  data.  This  matrix  should  be  used  by  those  who  want  to 
produce  a  rotation  that  is  as  close  as  possible  to  the  positions  determined  by  the  IERS. 
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Table  3.  Offset  angles  for  each  of  the  rotation  matrices. 


Matrix  parameters 
Equinox  ex,  ey  yy 

A  0 

Name 

ex 

(mas) 

ey 

(mas) 

yy 

(mas) 

A  0 

(mas) 

1 .  Inertial 

VLBI 

LLR 

LLR 

Riers 

-6.819 

16.6171 

40.83 

55.42 

2.  Inertial 

VLBI 

VLBI 

LLR 

-6.819 

16.6171 

38.328 

55.42 

3.  Rotating 

VLBI 

LLR 

LLR 

-6.819 

16.6171 

40.83 

-38.24 

4.  Rotating 

VLBI 

VLBI 

LLR 

-6.819 

16.6171 

38.328 

-38.24 

5.  Inertial 

LLR 

LLR 

LLR 

Rllr 

-5.36 

17.7 

40.83 

55.42 

Alternatively,  the  separation,  yy,  between  the  projection  of  the  equinox  onto  the  equator  of  the  ICRS  and  the  node  of  the 
ecliptic  on  the  equator  of  the  ICRS  may  be  determined  from  the  VLBI  offset  of  the  pole.  This  value  can  then  be  combined  with 
the  offset  between  the  equinox  on  the  ICRS  and  the  origin  of  the  ICRS  found  by  Chapront  et  al.  from  LLR  observations.  This 
change  in  the  derivation  of  yy  only  affects  the  R\2  and  Ro\  elements  of  the  rotation  matrix  as  given  in  (9). 

Similarly,  using  the  rotating  definition  of  the  mean  ecliptic  will  also  change  the  position  of  the  mean  equinox  relative  to 
the  origin  of  the  ICRS.  As  with  the  determination  of  yy,  the  use  of  the  rotating  definition  of  the  equinox  only  affects  the  R 12 
and  R21  elements  of  the  rotation  matrix.  Equation  (11)  gives  the  values  of  these  elements  using  the  rotating  definition  of  the 
equinox  with  ytJ  determined  from  LLR  observations,  while  (12)  gives  the  value  for  these  two  elements  with  yu  determined  from 
VLBI  observations. 

Finally,  the  matrix  /?llr,  (13),  is  based  solely  upon  the  parameters  determined  solely  from  LLR  data  by  Chapront  et  al.  This 
matrix  should  be  used  by  those  who  desire  a  rotation  based  upon  an  entirely  self-consistent  set  of  parameters  using  the  inertial 
definition  of  the  mean  equinox. 

Users  needing  accuracy  better  than  a  few  tens  of  mas  will  need  to  use  one  of  these  matrices  for  converting  from  J2000.0  to 
the  ICRS  and  the  subsequent  processes  of  applying  precession,  nutation,  and  the  formation  of  an  apparent  place.  However,  the 
user  needs  to  realize  that  there  are  historically  two  different  definitions  of  the  mean  equinox  that  differ  by  93.66  mas.  Finally,  the 
difference  in  the  VLBI  and  LLR  solutions  for  the  position  of  the  mean  pole  of  J2000.0  implies  that  there  may  be  an  inconsistency 
in  determining  the  equinox  that  can  lead  to  a  misalignment  in  the  origin  of  up  to  3  mas.  Thus,  care  must  be  used  in  selecting  the 
proper  matrix  to  rotate  between  the  ICRS  and  the  mean  dynamical  equator  and  equinox  of  J2000.0. 
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